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Abstract 

Q_^[ It appears difficult to construct a simple model for an open universe based on the 

one bubble inflationary scenario. The reason is that one needs a large mass to avoid 
the tunneling via the Hawking Moss solution and a small mass for successful slow- 
rolling. However, Rubakov and Sibiryakov suggest that the Hawking Moss solution 
rS I is not a solution for the false vacuum decay process because it does not satisfy the 



c^ I boundary condition. Hence, we have reconsidered the arguments for the defect of the 

simple polynomial model. We point out the possibility that one of the valley bounce 
belonging to a valley line in the functional space represents the decay process instead 
of the Hawking Moss solution. Under this presumption, we show an open inflation 
model can be constructed within the polynomial form of the potential so that the 
fluctuations can be reconciled with the observations. 
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1 Introduction 



Recent observations suggest the matter density of the universe is less than the critical density. 
Hence, it is desirable to have a model for an open universe, say Qq ~ 0.3. The realization of 
an open universe is difficult in the ordinary inflationary scenario. This is because if the universe 
expands enough to solve the horizon problem, the universe becomes almost flat. One attempt 
to realize an open universe in the inflationary scenario is to consider inside the bubble created 
by the false vacuum decay [|^. The scenario is as follows. Consider the potential which has 
two minimum. One is the false vacuum which has non-zero energy and the other is the true 
vacuum. Initially the fleld is trapped at the false vacuum. Due to the potential energy, universe 
expands exponentially and the large fraction of the universe becomes homogeneous. As the false 
vacuum is unstable, it decays and creates the bubble of the true vacuum. If the decay process is 
well suppressed, the interior of the bubble is still homogeneous. The decay is described by the 
0(4) symmetric conflguration in the Euclidean spacetime. Then, analytical continuation of this 
conflguration to the Lorentzian spacetime describes the evolution of the bubble which looks from 
the inside like an open universe. Unfortunately, since the bubble radius cannot be greater than 
the Hubble radius, the created universe is curvature dominated even if the whole energy of the 
false vacuum is converted to the energy of the matter inside the bubble 0. Thus, the second 
inflation in the bubble is needed. If this second inflation stopped when Q < 1, our universe 
becomes homogeneous open universe. 

Though the basic idea is simple, the realization of the scenario in a simple model has been 
recognized difficult 0. The difficulty is usually explained as follows. Consider the model involv- 
ing one scalar fleld. For the polynomial form of the potential like V{(j)) = m'^(f)'^ — 6(f)^ + A0^, 
the tunneling should occur at sufficiently large to ensure that the second inflation gives the 
appropriate density parameter. Then, the curvature around the barrier which separates the false 
and the true vacuum is small compared with the Hubble scale which is determined by the energy 
of the false vacuum. The fleld jumps up onto the top of the barrier due to the quantum diffusion. 
When the fleld begins to roll down from the top of the barrier, large fluctuations are formed due 
to the quantum diffusion at the top of the barrier. Then the whole scenario fails. The problem is 
rather generic. To avoid jumping up, the curvature around the barrier should be large compared 
with the Hubble scale V" > H^. On the other hand, to realize the second inflation, the fleld 
should roll down slowly, then we need V" < H^. These two conditions are incompatible. 

There are several attempts to overcome the problem. Recently Linde constructs the potential 
which has sharp peak near the false vacuum Q]. In this potential, the tunneling occurs and at the 
same time slow-rolling is allowed after the tunneling, then the second inflation can be realized. 
But, it is still unclear what is the physical mechanism for the appearance of the sharp peak in 
the potential. 

A more detailed study of the tunneling process is needed to tackle the problem. In the 
imaginary-time path-integral formalism, the tunneling is described by the solution of the Eu- 
clidean fleld equation. The solution gives the saddle-point of the path-integral. Then the solution 
determines the semi-classical exponent of the decay rate exp(— S'£;(0b)), where Se is the Euclidean 



action. In the case the curvature around the barrier is small compared with the Hubble, the so- 
lution is given by the Hawking Moss (HM) solution, which stays at the top of the barrier through 
the whole Euclidean time 0. Recently Rubakov and Sibiryakov give the interpretation of the HM 
solution in the de Sitter spacetime using the constrained instanton method |^, . They show the 
HM solution does not represent the false vacuum decay. This is because the HM solution does not 
satisfy the boundary condition that the field exists in the false vacuum at the infinite past. One 
should consider a family of the almost saddle-point configurations instead of the true solution 
of the Euclidean field equation. They show although the decay rate is determined by the HM 
solution, the structure of the field after tunneling is determined by the other configuration which 
is one of the almost saddle-point solutions. In the constrained instanton method, one introduces 
the constraint to define the subspace in the functional space and looks for the almost saddle-point 
solution. One must choose the constraint so that the almost saddle-point solution is contained 
in the region which is expected to dominate the path-integral. One way is to consider the valley 
region of the functional space |^, Q. Along the valley line, the action varies most gently. Then it 
is reasonable to take a configuration on the valley line as the almost saddle-point configuration. 
We will call the configuration on the valley line the valley bounce 0y. 

Their analysis opens up a new possibility to overcome the problem. Suppose that one of 
the valley bounces describes the tunneling and the inside of the bubble created from this valley 
bounce is our open universe. The initial condition of the tunneling field in the open universe is 
determined by the valley bounce. If the field appears sufficiently far from the top at the nucleation, 
the large fluctuations can be avoided. During the tunneling, fluctuations of the tunneling field 
are generated. These fiuctuations are stretched during the second infiation and observed in the 
open universe. 

In this paper, we extend the valley method developed by Aoyama.et.al to the de Sitter 
spacetime. Taking the assumption that one of the valley bounces describes the tunneling and 
inside the bubble created from the valley bounce is our universe, the fiuctuations generated during 
the tunneling are calculated by defining the fiuctuations which are relevant to the observable in 
the open universe as those orthogonal to the gradient of the action. Then we show an open 
universe can be constructed from the valley bounce and the fiuctuations can be generated with 
the appropriate properties. The paper is organized as follows. In the next section we review 
the formalism to describe the false vacuum decay in the de Sitter spacetime. Then we explain 
the role of the configurations on the valley of the action, i.e. valley bounces. We derive the 
valley equation which determines the valley bounce. In section 3, we solve the valley equation 
analytically using the piece-wise quadratic potential and fixed background approximation. Then 
the structure of the valley is shown. In section 4, we explain the scenario emphasizing the role of 
the valley bounce and show that a simple model for an open universe can be constructed without 
introducing the fine-tuning of the potential. In section 5, we calculate observable fiuctuations in 
the open universe from the valley bounce. The power spectrum of the curvature perturbations is 
calculated and these fiuctuations are shown to be compatible with the observations. In section 
6, we summarize the results. 



2 Valley method in de Sitter spacetime 



First we review the formalisms which are necessary to describe the false vacuum decay in 
the de Sitter space. We want to examine the case in which the gravity comes to play a role. 
Unfortunately, we have not known how to deal with quantum gravity effect yet. So, we study 
the case in which we can treat gravity at the semi-classical level. That is, we treat the problem 
within the framework of the field theory in a fixed curved spacetime . The potential relevant 
to the tunneling is given by 

V(0) = e + VV(0). (1) 

We assume e is of the order Mf and Vr(0) is of the order M'^. We study the case M is small 
compared to M^, M -C M^. Then the geometry of the spacetime is fixed to the de Sitter spacetime 
with H = M^/Mp , where M~^ = SnG/S. We consider the situation in which the potential Vt{4>) 
has the false vacuum at = ifp and the top of the barrier at (p = ifT- Since the background 
metric is fixed, we can change the origin of the energy freely. We choose Vr(0F) = 0. Following, 
we work in units with H = 1. 

The decay rate is given by the imaginary part of the path-integral 

Z = J[d4>]exp{-SEm, (2) 

where Se is the Euclidean action relevant to the tunneling. The dominant contribution of this 
path- integral is given by the configurations which have 0(4) symmetry [^. So, we assume the 
background metric and the field to have the form 

ds'^ = da"^ + a(o")^ (dp^ + sin^ pdQ,) , 
cP = 0(a), (3) 

where a{a) = sin a. Then, the Euclidean action of 0(cr) is given by 

Se = 27r2 J da (a' (^<P'' + Vt{<P)^ ) . (4) 

The saddle-point of this path-integral is determined by the Euclidean field equation 6Se/S(I)=0; 

0" + 3cot(T(/)'-y^(0) =0. (5) 

We impose the regularity conditions at the time when a{a) = as 

0'(ct = 0) = (p'{a = tt) = 0. (6) 

We represent the solution of this equation as (psi.'^)- If the fluctuations around the solution have 
a negative mode, it gives the imaginary part to the path-integral and this solution contributes to 
the decay dominantly. The decay rate F is evaluated by 

F~exp(-^s(0B)). (7) 



The equation has two types of the solutions depending on the shape of the potential. If the 
curvature around the barrier is large compared with the Hubble scale, then the Coleman De 
Luccia (CD) solution and the Hawking Moss (HM) solution exist p, [1^. Since the CD solution 
has lower action than that of the HM solution, the decay is described by the CD solution. The 
analytic continuation of this solution to Lorentzian spacetime describes the bubble of the true 
vacuum. On the other hand, in the case the curvature around the barrier is small compared with 
the Hubble scale, only the HM solution exists. This solution is a trivial solution cj) = ipx- The 
meaning of the HM solution is somewhat ambiguous. There are several attempts to interpret this 



tunneling mode. One way is to use the stochastic approach ^^. It has been demonstrated that 
the decay rate given by eq.(7) coincides with the probability of jumping from the false vacuum 
i^F onto the top of the barrier ipx due to the quantum fluctuations. 

Recently, Rubakov and Sibiryakov gave the interpretation of the HM solution using the con- 
strained instanton method [§, . The main idea is to consider a family of the almost saddle-point 
configurations instead of the true solution of the Euclidean field equation, i.e. the HM solution. 
The motivation comes from the boundary condition. They take the boundary condition that 
the state of the quantum fluctuations above the classical false vacuum is the conformal vacuum. 
For this boundary condition they show the field should not be constant at < a < vr and the 
HM solution is excluded by this boundary condition. Then one should seek the other configura- 
tions which obey the boundary condition and dominantly contribute to the path-integral. In the 
functional integral, the saddle-point solution gives the most dominant contribution, but the con- 
tribution from a family of almost saddle-point configurations which have almost the same action 
with that of the saddle-point solution should also be included. To seek the almost saddle-point 
solution, one introduces the constraint to the path integral. The constraint selects the subspace 
of the functional space. The minimum in this subspace satisfies the equation of motion with 
constraint instead of the field equation. This minimum corresponds to the almost saddle-point 
configuration which is slightly deformed from the HM solution. Since the HM solution gives the 
minimum action, the decay rate is determined by the HM solution. But the structure of the field 
after tunneling can be determined by the one of the almost saddle-point configuration. They 
found that the configuration describes the bubble of the true vacuum in the Lorentzian space- 
time. Then, they conclude that even in the case only the HM solution exists, the result of the 
tunneling process can be the bubble of the true vacuum which is described by one of the almost 
saddle-point configurations. 

In the constrained instanton method, the validity of the method depends on the choice of the 
constraint O, O]. We should choose the constraint so that the almost saddle-point solution is 



contained in the region which is expected to dominate the path-integral. Since the action varies 
most gently along the valley line, we consider the valley region of the action |^, Q. 

In the false vacuum decay process in the fiat spacetime, the configurations along the valley 
line have physical meanings [|14| . These configurations actually dominate the path-integral in the 
presence of the low energy incoming particles. If the incoming particles exist, the path-integral 
which describes the decay is given by /[(i0]0'^exp(—S'E ((/>)). The effect of the incoming particles 
deforms the saddle-point. As long as the energy of these incoming particles is low, the deformed 
configurations belong to the valley. This is because in deforming the configurations from the 
saddle-point solution, the configurations on the valley line can be obtained most easily compared 



with the other configurations with the same action. Thus, the configurations on the valley play 
a crucial role to calculate the decay rate or the cross section when the initial state is of higher 
energy than the ground state. In the de Sitter space time, the choice of the quantum states of 
the quantum fiuctuations above a given classical false vacuum will affect the tunneling process 
considerably |^. Hence, we think the configurations on the valley line in the de Sitter spacetime 
may also play an important role, though it is not easy to identify the corresponding quantum 
state. 

Taking into account the above fact, it is desirable to analyze the structure of not only the 
solution of the Euclidean field equation but also the configurations on the valley line. One way to 
fined the configurations on this valley line is to use the valley method developed by Aoyama.et.al 
1^. To obtain the intuitive understanding of this method, consider the system of the field 0j. 
Here i stands for the discretized coordinate label and we take the metric as Sij. In the valley 
method the equation which identifies the valley line in the functional space is given by 

D^jdiS = XdiS, Dij=d^djS, (8) 

where di = d/d(j)i. Since the equation (8) has one parameter A, the equation defines a trajectory 
in the space of 0. The parameter A is one of the eigen value of the matrix Dij. On the trajectory 
the gradient vector diS is orthogonal to all the eigenvectors of Dij except for the eigenvector of 
the eigen value A. The equation can be rewritten as 

5.Q(5.5)'-A5)=0. (9) 

Then the solution extremizes the norm of the gradient vector diS under the constraint S =const., 
where A plays the role of the Lagrange multiplier. Such solution can be found each hypersurface 
of constant action, then the solutions of the equation form a line in the functional space. If we 
take A as the one with the smallest value, then the gradient vector is minimized and the action 
varies most gently along this line. This is a plausible definition of the valley line. We will call the 
configuration on the valley line of the action the valley bounce (pv and the trajectory they form 
the valley trajectory. 

We shall formulate the valley method in the de Sitter spacetime. The most convenient way is 
to use the variational method eq.(9). We shall define the valley action by 

'^-'^-lxl''^{y-!f)^- (10) 

The valley bounce is obtained by varying the action Sy- The equation which determines the 
valley bounce SSy/Scp = is a fourth order differential equation. We introduce the auxiliary field 



/ to cancel the fourth derivative term [|14 



''-rJ''^^^(^-jlfl (" 



Then the valley action becomes 

2XJ ^^' A 



1 /■ 1 /■ (^9 

Sv + Sf = SE + ^ da^f - - / daf-§. (12) 



Taking the variation of this action with respect to / and 0, we obtain the equations for and /; 

1 6Se 



v^ 



/, 



6^Sf 



Using a{a) = sin cr, the valley equation which determines the structure of the valley bounce is 
given by 

0" + 3 cot a 0'- 14(0) = -/, 
r + 3cota/'-l4'(0)/ = -A/. (14) 

The fluctuations around the valley bounce can be expanded by the eigenmodes ga,n{o') with the 
eigenvaluepa,n of the operator {6S'^/6(f)6(t>)^^; 

9'a,n + 3 cot a g'^^^ - V^{(j))ga,n = -Pa,nga,n- (15) 

Since A is one of the pa „ with the smallest value, the gradient of the action f{a) is orthogonal to 
the other eigenmodes with p„ „ 7^ A. To ensure that the valley bounce gives the imaginary part to 
the path-integral, the fluctuations around the valley bounce should have one negative eigenvalue 
P- 

In the next section, we solve the valley equation and clarify the structure of the valley bounce. 
Solving the valley equation is the eigenvalue problem of the two variables, it is desirable to solve 
the equation analytically to confirm the existence of the solutions. The valley bounce can have 
a thick-wall profile. In the fiat spacetime, there exists the attempt to treat thick-wall solutions 
analytically by constructing the piece- wise quadratic potentials |jlll. In the next section, we 



extend the attempt to the de Sitter spacetime and solve not only the Euclidean field equation 
but also the valley equation analytically. 



3 Valley bounces 

3.1 The construction of Valley bounces 

To solve the valley equation developed in the last section, we construct the piece-wise quadratic 
potential. We connect two parabola. In the potential the true vacuum is absent. But this is 
not essential in calculations. In fact we have solved numerically the valley equation for several 
potentials and found this model is sufficient to discuss the generic feature of the valley bounce. 

The potential which we study is 

^m%{(f)- (PfY, -00 < < 0, 

Vt{<P) = { (16) 

-\m\{(t) - (^tT + r], < < 00, 

7 



where 77 is of the order M^. We require that the potential and its derivative are connected 
smoothly at the connection point = 0. From this condition, we obtain 



m. 



ifT 



Vf 



rrin 



-Vf, 



2my 1] 



ipT and ifF have a mass scale of the order hP Jvht- Thus, we rescale the field as 
and / -^ {M^/mT)f. 

First we solve the Euclidean field equation; 

0" + 3 cot a 0'- 14(0) =0. 



(17) 



[li 



If one puts z = —cosa and Y{z) = Vl — -2^(0 — V^i), the equation reduces to the associated 
Legendre differential equation 



'l-z' 



(PY 



dV 
dz"^ dz 



Vi{Vi + 1) 



^^ 



1-^2 



r = o. 



where n = 1 and Ui is given by 



/9 



1 



ut 



- + m2.--. Up 



-mp - 



(19) 



(20) 



Here, i = T for —00 < < and i = F for < < 00. The independent solutions of the 
equation are given by the associated Legendere function of the first and second kinds, Pj, {z) and 
Ql.{z). Pl{z) is regular at 2; ^ 1. Since these functions behave at 2; — > — 1 as 



the combination of these solutions 



-2^/'sin(7rz/)7r-^(l + 2)-^/^ 
-2-i/2^os(7rz/)(H-2)-^/^ 



B^^{z) = P!:iz) + (^--tan(7rz/)j Q^^iz) 

is regular at z ^ —1. Then the solution which satisfies the boundary condition is given by 

1 



(21) 



(22) 



>i?= < 



(fF 



ipF + 



1 



AFBl^iz), 



AtPIM^ 



-1 < Z < Zq, 



Zq < Z < 1, 



(23) 



where 0b (-Zq) = 0. Since the potential is constructed to be smooth to its first derivative, we 
demand 0^ and its first derivative must be continuous aX z = zq. Then the coefficients Ai are 
determined in terms of zq; 



Af{zo) 



JfiF V 1 - ^0 



At{zo) 



ifT^l-Z^ 



(24) 



The junction time zq is determined by 0b(-2o) = 0; 



(25) 



If the algebraic equation for zq has a solution, this gives the CD solution. The condition for the 
existence of the solution restricts the parameter rrii. We see this condition is approximately given 
by m| > 4. 



Next we shall solve the valley equation. Equation for / is given by 

r + 3cota/'-\4'(0)/ = -A/. 



(26) 



The regularity conditions are the same with that of 0. Then, the general solution which satisfies 
the boundary conditions is given by 



/ = 



VT 



GpBl (z 



I VT^^ 



upx V / ' 



GtPI,,{z) 



where 



I'TX 



'9 1 



-l<z< Zx, 
z\< z<l, 



\-{ml-\)-\. 



(27) 



(^Ti G F 



PlAzx) 



l^TX 



Cjt- 



From the junction conditions, we obtain 

'- Bl^^{z,y - B^^iz,) 
The equation has solutions only if zx satisfies the following equation 

Pl,Szx)Bl^izx) - Pl^izx)Bl^izx) = 0, 
which determines the junction time zx- Next we solve the equation for 0; 

0" + 3cota0'-\4(0) = -/. 
In the equation, / acts as the source. We can see that the special solution is given by 

/ 



(p-Vi 



X' 



f28) 



(29) 



(30) 



(31) 



(32) 



Then the solution which satisfies the boundary condition is given by 



v^ 



AVI 



^T + ^=i= At PUz) + -^=L= Gt PLAz), 



VT^ 



XVT^Y^ 



-1 < Z < Zx, 

zx<z<l. 



(33) 



From the junction conditions, we obtain the coefficients 

'^' ''PL(zx)Blh,)-PL(zx)Bl(z,] 



Ap = \l~z 



l/^ V '^ / Z^i? V '^ / V'T^'^-'i^F^ 



G^ = -\,ll-z\ VfPIM)BIA^>^)-VtPU^x)BI^{z,) 



^ Bl^^{z,){Pl^{z,)Bl^{zy) - Pi,{z>)Bl^{z>)y 

Note that in this model the deformation of the configurations is essentially determined by z\. If 
zx = zo, Gi becomes 0, so cj)v = (pB as expected. 

To ensure that the valley bounce plays a role instead of the true saddle point solution, we 
must examine the fluctuations around the valley bounce have one negative mode and give the 
imaginary part to the path- integral. The equation which determines the eigenmodes is given by 

g'^ + 3cota g'^- VJ^{(t))gn = -pngn- (35) 

Then, the eigenvalue equation which determines the eigenvalue pn of these eigenmodes becomes 

Plr.^ i^x)Bl,,^ (^a) - Pl^^ {zx)Bl^^ {zx) = 0. (36) 

Note that, A is one of the solutions p„. 

We should treat separately the case in which the valley bounce exists around the top of the 
barrier and passes through only one parabola. We put the solution for / as 



oo 



f = J2^nCosna, (37) 

n=0 



then the equation for / is rewritten as [T^ 



J2 {[{n - l)(n + 2) - m| - A]6„_i - [{n + l)(ra - 2) - m^ - A]6„+i) sinna = 0. (38) 



n=0 

Thus, bn converges only when 

X = -ml + n{n + 3). (39) 

We put the solution for as 

oo 

— (/9t = ^ a„ COS na, (40) 

n=0 

then the solution for is given by 

On = jCn- (41) 

The eigenvalue of the eigenmode at the valley bounce is given by 

p„ = -m| + n(n + 3). (42) 
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3.2 The structure of the valley 

Using the analytic solution of and /, we show the structure of the valley. Remember that we 
have rescaled the field as — >■ {M^/mT)(f) and / — *> {M"^ /mTjf ■ Following, for completeness, we 
consider the two types of the potential; {l)m\. > 4 and (2)m^ < 4. 

(l)mf, > 4 

There exist two solutions in the Euclidean field equation; the CD solution and the HM solution. 
For example, we take m^ = 7, mjp = 2.2 and r] = 0.6M^ (Fig.l). The behaviors of the CD solution 
0(o"), the eigenmode with the negative eigenvalue g-{(y) and the scale factor a{a) are shown in 
Fig. 2. The CD solution has one negative eigenvalue pco.- = —4.7 and the smallest positive 
eigenvalue is pcD,+ = 3. Since the CD solution gives the saddle-point of the path-integral, we 
analyze the valley trajectory which contains the CD solution. At the CD solution on the valley 
trajectory, / = 0. The valley bounce near the CD solution is obtained by deforming the CD 
solution; (py = (pcD + A0y. The deformation Acpy is due to the source term / 7^ in the 
equation for 0. The equation for / is almost the same with that for the eigenmode g at the CD 
solution. Thus, A is given by A = pcD + AA. To ensure the action varies most gently along the 
valley trajectory, pcD should be the eigenvalue of the smallest value. In this case, pcn has one 
negative eigenvalue, so we take A at the CD solution as X{4>cd) = Pcd- or A(0(7d) = Pcd,+- 

First examine the valley trajectory associated with the negative eigenvalue (A(0(7d) = Pcd-)- 
The valley bounce obtained from the analytic results developed in the previous section is shown 
in the lower-panel of Fig. 3. Since pen- is the lowest eigenvalue, / does not have a node. In 
the equation for 0, / acts as the force. So, the valley bounce in this trajectory is obtained by 
deforming the CD solution adding a one-direction force /. If / > 0, the valley bounce has a 



structure of the small bubble and if / < it has a structure of the large bubble |jTj] . We plot the 
action along this trajectory in lower-panel of Fig.4. The CD solution gives the maximum of the 
action. 

Next consider the valley trajectory associated with the smallest positive eigenvalue (A(0cd) = 
pcD,+)- We show the valley bounce in this trajectory in the upper panel of Fig. 3. Since pcD,+ 
is the next to the lowest eigenvalue, / has one node. In the equation for 0, / acts as the mass 
term of 0. So, in this trajectory, the valley bounce is obtained by modifying the mass of the field 
0. It is known that the CD solution is smoothly connected to the HM solution if one decreases 
the mass around the top of the barrier [|T^. Then, it is expected that this trajectory connects 
the CD solution and the HM solution. The action along this valley trajectory is shown in the 
upper-panel of Fig.4. Since the degeneracy occurs in A, we take the horizontal coordinate as the 



'norm' of the solution |$| = J 271"^ J a^ {(f) — (PtY [jT4[- We see the CD solution is the minimum 
and the HM solution is the maximum of the action and these solutions are smoothly connected 
on this trajectory as expected. 

(2)m^ < 4. 

The saddle point solution is the HM solution. For example we take m"^ = 2, ml. = 0.5 
and rj = O.IM^ (Fig.5). The HM solution has one negative eigenvalue Phm- = —2 and the 
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smallest positive eigenvalue is given by Phm,+ = 2. The generic feature of the valley bounce 
is understood by the simple analysis of the case in which the valley bounce exists only in one 
parabola. First consider the valley trajectory associated with the negative eigenvalue. The 
solution of the valley equation is essentially has a form / = A(0 — ipt) = const. This solution 
does not represent the tunneling, so we seek the trajectory associated with the smallest positive 
eigenvalue X{(Phm) = Phm,+- The solution of the valley equation is given by — 0^ oc coscr 
and / = X{(j) — (pr) (Fig-6). In this trajectory, the HM solution gives the minimum of the action 
(Fig. 7). The action grows as the variation of the field becomes large, but the increase is relatively 
gentle. 

The fluctuations around the valley bounce should have one negative mode to ensure that the 
valley bounce plays a role instead of the HM solution. The valley bounce has a lowest eigenvalue 
Pv- < A(0v'), which is negative on this trajectory. Since this is the unique negative eigenvalue, 
the gaussian integration of the fluctuations around this valley bounce gives the imaginary part to 
the path-integral. Then, the valley bounce contributes to the false vacuum decay and describes 
the creation of the bubble of the true vacuum. 



4 An open universe from valley bounce 

Using the results developed so far, we will study a model for an open universe inside the bubble. 
Since the radius of the bubble R is small compared with the Hubble horizon |]TU[, then the 
curvature scale is greater than the energy of the matter inside the bubble pm even if the whole 
energy of the false vacuum is converted to the energy of the matter, pm/M^ ~ H^ < 1/ E? 0. 
Thus, we need the second inflation in the bubble. To realize the second inflation inside the bubble, 
the field should roll slowly down the potential. It requires that the curvature of the potential 
is small compared with the Hubble. To avoid the ad hoc fine-tuning of the potential, we will 
assume the requirement is satisfied for all region of the potential. Since ttit < H, the solution 
of the Euclidean field equation is given by the HM solution. If the tunneling is described by the 
HM solution, the field appears at the top of the barrier. Then large fluctuations are generated 
because at the top of the barrier the field experiences the quantum diffusion rather than the 
classical potential force. Fluctuations in this diffusion dominated epoch make the inhomogeneous 
delay of the start of the classical motion, thus make large fluctuations. 

The above argument is based on the assumption that the HM solution describes the false 
vacuum decay. However, it seems to be possible that one of the valley bounces describes the 
decay instead of the HM solution. One of the grounds is that in the de Sitter spacetime, the 
choice of the quantum states of the quantum fluctuations above the false vacuum will affect the 
tunneling process considerably 0. In the de Sitter spacetime, one has to specify the state of 
the quantum fluctuations of the field besides specifying the classical vacuum which is the average 
value of the field i.e. (fp- The dominant configuration in the path-integral depends on the initial 
state of the quantum fluctuations above the false vacuum. In case of the flat spacetime, if the 
initial state is of higher energy than the ground state, the dominant configuration is given by one 



of the valley bounces instead of the bounce solution [14|. Thus it seems natural to consider the 



situation in which one of the valley bounces describes the tunneling and inside the created bubble 
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is our universe, although it is difficult to identify the quantum state corresponding to the valley 
bounce. Hence we take the presumption that one of the valley bounces describes the tunneling 
and inside the bubble created from the valley bounce is our universe. 

A problem about the assumption that an individual valley bounce describes the tunneling is 
the interpretation of the valley bounce in the Lorentzian spacetime. Our interpretation is that the 
valley bounce determines the initial condition of the tunneling field in the open universe. Within 
the fixed background approximation, we can make analytic continuation about the background 
geometry from the Euclidean de Sitter space (3) (see Fig.8). By the analytic continuation 

r = 2(p-7r/2), a = a, (43) 

we obtain the Lorenzian de Sitter spacetime (Region II) 

ds"^ = da^ + a{af (-dr'^ + cosh^ rd^^) . (44) 

We take the nucleation surface at r = 0. Region II is almost covered by the false vacuum. Then 
we assume the effect which modifies the dominant configuration from the HM solution to the 
valley bounce also modifies the classical motion of the field in this region. The field obeys the 
equation analytically continued from the valley equation (14). The solution of the equation is 
given by the analytic continuation of the valley bounce. On the other hand, because the bubble 
expands classically at a velocity rapidly approaching the velocity of light, inside the expanding 
bubble is well described by the usual classical equation of motion. On the light-cone of the center 
of the bubble a = vr, the coordinate is singular a{a) = 0. We continue to the interior of the 
light-cone (Region I) by 

r = T + i—, t = i{a — n). (45) 

The resulting metric is given by 

ds^ = -dt^ + h{tf (dr'^ + sinh^ rd^^) , (46) 

The expanding bubble is homogeneous and isotropic on the hypersurface on the hyperbolic time 
slicing t = const. The interior of the light-cone can be viewed as an open Friedman-Robertson- 
Walker universe with scale factor b(t). The initial condition of the tunneling field on t = 
hypersurface is determined by the behavior of valley bounce at c = vr. After that time, the 
evolution of the field is described by the classical field equation. 

Under these assumptions, we will construct a model for an open inflation in the simple model 
with the polynomial form of the potential. We connect the linear potential at the point the field 
appears after the tunneling = 0*, 

l^(0) = K-/i'(0-0*), (0>0*). (47) 

We demand the potential and its derivative are connected smoothly at the connection point 0*. 
Then we obtain 

K = e + ?7 - -m|(0* - v9t)^, 

/i^ = m|(0, -v?r)- (48) 
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The initial conditions of the field in the open universe are given by the valley bounce 

0(t = 0) = 0o(^ = 1) = 0*, 0(t = O)=O. (49) 

The field evolves obeying the classical field equation; 

+ 3cotht0 + V^'(0) = 0, (50) 

then the solution of (h satisfies 



3 cosh t — 3 cosh t + 2 



^W^" siiS?^ ^ <^i' 

In the small t this behaves as (l/4)/x^t. The classical motion during one expansion time is given 
by |0|iJ^^. On the other hand the amplitude of the quantum fiuctuations is given by 50 ~ H . 
The curvature perturbation TZ produced by the quantum fluctuations is approximately given by 
the ratio of these two quantities; 

5<p H^ H"^ ( H \ 

n ~ ^— ^ . 52 

This should be of the order 10^^ from the observation of the cosmic microwave background (CMB) 
anisotropies. If |0* — ^^rl < -f^, as in the case the HM solution describes the tunneling, IZ > 1 
and the scenario cannot work well. Fortunately, from Fig. 7, we see for appropriate A, the valley 
bounce gives the initial condition as |0* — ipx\ ~ 0(l)(M^/mT'), which is larger than the Hubble 
ii M > H. For this initial condition, the potential force works and the fleld rolls slowly down the 
potential according to eq.(50). We expect the curvature perturbation can be suppressed for the 
valley bounce 



5 Fluctuations in the open universe 

In this section, we will calculate the observable fluctuations in the open universe from the valley 
bounce based on a model described in the previous section. We will assume the decay is described 
by one of the valley bounces 0o- Then we take one speciflc A in the following calculations. Under 
this assumption, we calculate the observable fluctuations in the open universe. 



5.1 Fluctuations around the valley bounce 

We flrst calculate the fluctuations around the valley bounce in the Euclidean spacetime. The 
Euclidean action can be expanded around the valley bounce 0o as 



5'e(0) = S'i5(0o) + j d^xy/g— -jj- 



^) + l f fd^^d'x' ^'^' 



2J J 6(j){x)6(f){x') 



6(l){x)6(l){x'). (53) 
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Since the valley bounce does not obey the field equation, then the first order derivative term does 
not vanish. This can be avoided by constraining the space of the fiuctuations to that orthogonal 
to the gradient of the action; 



,4 ^ / 1 SSf 



.V9 



(54) 



We must calculate the physical observable like two-point correlation function in the bubble 
described by 0o- For example, consider the variance of the scalar field (0^) — (0)^, where ( ) is 
average over p and Q. Analytically continuing to the Lorentzian spacetime, this corresponds to 
the average over the space in open universe. Since the variables which depend only on a obey 
the relation 

(0o)=0o, (<50(a)) = 50(a), (55) 

we can show 

(02) - (0)2 = (50(a, p, m - (<50(a, p, n)f. (56) 

So the observable in the open universe can be evaluated from inhomogeneous fiuctuations which do 
not have 0(4) symmetric configurations. In the de Sitter spacetime, inhomogeneous fluctuations 
are expanded by scalar harmonics 

6(j){a,p,Q) = I dpSp{cr)Ypi^{p,Q). (57) 

The harmonics obeys the orthogonal relation between different p^. Now the gradient vector is 
given by /(cr) which is the mode of p^ = — i. Then, the inhomogeneous fluctuations which depend 
on p and Q {p'^ ^ —1) are orthogonal to the gradient of the action f{a) automatically. 

Fluctuations around the valley bounce 0o obey the field equation 



dx''^ 



6'^Sf 



x') = 0. (58) 



6(f){x)6(f){x') 

In Region II of the Lorenzian de Sitter spacetime; 

ds'^ = da'^ + a{af (-dr"^ + cosh^ rdn^) , (59) 

we assume background field obeys the equation of motion analytically coninued from the valley 
equation. Then the fiuctuations in Region II is obtained by the analytical continuation (45); 

r = z(p-7r/2), (j = a. (60) 

Thus we will solve the equation (58) analytically continuing to Region II. The procedure to solve 
the equation is the same with that was done in the previous works [T^, |TH|, [T^, Q . We will follow 
their calculations. Expanding the fiuctuations as 



50(cr,r,fi) = / dpSp{a)Ypi^{T,VL), (61) 
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we obtain the equation of the fluctuations 



where 



f — + 2tanhr— + -^-^jy;,„(r,f]) = -(l+/)y;,„(r,fi), (62) 

S;{a) + 3cotaS';(a) + (^^ - \4'(0o)l S,{a) = 0, (63) 

{ml, < cr < (Ta, 

— -m^, O'A ^ o" ^ TT, 

and 2;^ = — coso-^. Here A is determined by 0o- Since the temporal coordinate r is included in 
the harmonics l^/m, the choice of the solution Ypi^ specifies the vacuum. This choice is related 
to the initial quantum state of the fluctuations. We will take this initial state as Bunch-Davis 
vacuum. The equation for Sp{(j) can be rewritten as 

where 

a{u) = (cosh-u)~^, U{u) = j 5 tanh-u = — coso" = z. (65) 

cosh u 

Since U{u) ^ as -u ^ ±oo, the modes are continuous for p^ > 0. For ux < u, the potential has 

a valley, then some discrete modes exist for p"^ < 0. 

First consider the continuous modes. Positive frequency mode should satisfy the Klein-Gordon 
normalization 



— i dz 



COsh^ r dQ (50+,„ (<9r<50+,;,„,) - (5r'^0pZ„) S(pp*l'm')] _„ = S{P - P%l'Smm'- (66) 



T = 

For simplicity we consider s-wave. The normalized positive frequency mode function of the 
Bunch-Davis vacuum is given by 

■Kp/2 -ipT _ -■Kp/2 ipT 

5<Plp{a, r) = S^pia)QpiT), Qpir) = _^==^-___, (67) 



where Sp is normalized as 



' dz Spiz)S;,iz) = -^5ip-f). (6J 



-1 ^^~'-p'^~' SttH 

Using this mode function, the fluctuations can be expanded as 

5(j)= dp {5(j)'^ CLp + 5(j)tp a^p) + {h.c) , (69) 

•^ u 

where cip annihilates the Bunch-Davis vacuum. We take the initial fluctuations Sp{z) at z = — 1 
as the Klein-Goldon normalized mode Fp[z) on — 1 < z < 1, then we evolve this mode using the 
field equation to 2; = 1. The resulting mode function is 

if W, -l<z< ZA, 

Sp{z) = { (70) 

apF;^{z) + ppFlpiz), zx<z<l, 
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where 



and 



1 1 

FHz) = - 



a\T{l - tp)Pt^{z) - air(l + tp)p- 



ip( 



(71) 



\ 



i + \/i-|C|IVQ' 



ci 

\a\ 



\ 



l-Jl-\C\\^/Cf 



Clip) = 2vr 1 + ^^ , C^,{p) = -2n 
\ smn Tip J 



2r[l — ip] sinvrz/j 



r[l + ip] sinh Tip T[—ip — z/i]r[l — ip + Vi 



Here, «„ and /3p are determined by the junction conditions at zx 



ap{zx) 
(3p{zx) 



J?Fj pT _ pT ^ pF 



^p (''zJ^ —p ^ —p^z^ p 

-F^d F^ + F'^d F^ 



J^ p ^Z^ —p ^ —p^z^ p 



2a 



z\ 



(72) 



where dz = d/dz. Sp{z) is not normahzed on — 1 < 2; < 1. The normahzed mode function is 
given by 



Sp{z) 



b+F^{z) - b^FLiz 



p\-j, ~l < Z < Zx, 



(73) 



(b+ap - b^P_p)F^{z) + {b+Pp - b_a_p)F^p{z), zx < z < 1, 



where 

b+ = 
and 



Di 



D'i 



I n |2 



\ 



1 + Jl - ID2PM 



6_ 



D2_\ / Di . 

ID2I; V^?-|^2|M 



1- J1-ID2IVA 



(74) 



Di{p) = -(|5p|2 + |/3p|2 + l), D2{p) = &p^p + 



cl_ 

2Cr 



Cir, 



KjLiqKXi \ kJifyiXi • r~^X) — r-^'D -\- LX'pLt _ 



(75) 



Next consider the discrete modes. We put p^ = — A^. The Bunch-Davis positive frequency 
mode is given by 

S(l)tirr. = SAia)YMrniT,n), (76) 

where 



-1-1/2, 



YMm{r, Q) 



and Sa is normahzed as 



^r[A + / + i]r[-A + / + 1] PA:i;r(^sinh 

2 Vicoshr 



dz\SA{z)\' = l. 



T 



(77) 
(78) 
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From the regularity condition similar to the valley bounce, the solution is given by 



^a(^) 



Vi^ 



{ptA^) + pKP;,\z)), 



[ Vl^z 



^P.-J'i^), 



-1< Z < Zx, 
zx< z <1, 



where 



/9a = ^^^^^^r[l + A + UF]r[A - up]. 



n 



From the junction condition, a\ is given by 



P. 



CtA 



ux 



Pt + ^IvP.-'' 



Up 



Up 



"A 



dzP,7^ 

Z l/rp 



Z\ 



d,Pl + /?A4Pr^ 



Up 



2a 



Then A is determined by the equation 

p,;^(4p4 + /5a4P-/) = d^P-J'iPup + PaPuJ") 

The mode Sa is not normalized. The normalized mode is given by 



Sa{z) = NaSa{z), Na 



1 . \ -1/2 

dz\SA{z)\' 
1 



(79) 



^0) 



f82) 



^3) 



5.2 Initial fluctuations in the open universe 



Fluctuations propagate into the interior of the light-cone a = 'k{z = 1). Since the coordinate 
system (59) is singular on the light-cone, we make analytic continuation by 



TT 



r = T + i—, t = i{a — TT). 

The resulting metric is given by 

ds^ = -dt^ + h{tf (dr^ + smh\dn^ 



U) 



^5) 



where b{t) = sinht. Since the fluctuations exponentially expand during the second inflation in 
the bubble, the shortwavelength modes are relevant. The matching condition across the lightcone 
in the Minkowski limit is given by 






2^J^ y/2 sinh np 
-i 1 



2^y2p v^2 sinh np 






m) 



where 



Uv) 



Rpir) 



-tprj—T] 



, e"^ = tanh(t/2). 



1 sin p r 
a/Stt sinh r 



^7) 



Note that Rp{r) is the norinahzed scalar harmonics Rp{r) = l^oo(^)) where 

YMr,^) = '^I^^^^^^^^^^YUni (88) 

and Yim is the usual spherical harmonics. Then, the extension to the general modes with / 7^ 
0,m 7^ is straightforwardly given by replacing Rp{r) to Ypimir^Vt). We obtain the fluctuations 
inside the bubble 

X [((e-P/^^i(p, A)T,(r^) + e--^''g,{p, X)T_,{r^))a, 

+ {e-^'^gliv. ^)T^M + e-'^gliv. A)T,(r7))a_p) +(/i.c.)] , (89) 

where 

5-1 (p, A) = a^^{h+ap{zy) - b_p^p{zx)) - a^*{b+pp{zx) - b^a^p{zx)), 

g2{p, A) = a1{b+Pp{zx) - b^a^p{zx)) - a^_{b+ap{zx) - b-(3-p{zx)). (90) 

This is the initial condition of the fluctuations in the open universe. The discrete mode can be 
treated in the same way. We obtain the positive frequency mode 

HI = Na ^^^^^^^^^YuUr, fi). (91) 

smht 



In the limit t ^ 0, this becomes 



rTAiri)YMm{r,^), (92) 



where Ta = e^^~^. Using this mode function, the fluctuations can be expanded as 



Y.Y.HlAA. + ih.c). (93) 



Im 



In Fig. 9 we plot the solutions A for m^ < 4 (Case (2) in section 3.1). We also show the normal- 
ization factor Na- We flnd two solutions of A. We call the mode with < Agub < 1 the subcritical 



mode and the one with 1 < A^^p the supercritical mode [^ . In the case the background solution 
is given by the HM solution, one supercritical mode with Asup = vt = \/9/4 + mf^ — 1/2 > 1 
exists. In the present case, the mass changes rn^ to —rn^ at zx-, another subcritical mode appears. 
Note that in the case the CD solution describes the tunneling, the supercritical mode corresponds 
to the wall fluctuation mode 50^ oc (f)cD with Kg^p = 2. Although in the present case the corre- 
spondence cannot be held, the behavior of the supercritical mode resembles the wall fluctuation 
mode. 
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5.3 Curvature perturbations in the open universe 

In this subsection we restore the Hubble scale H. The field evolves with the classical field equation 
inside the bubble (Region I). We should match the fiuctuations (89), (93) to the solution of the 
field equation with the background field satisfying the classical equation of the motion (50). 
Fluctuations of the scalar field give rise to a metric perturbations in the open universe. So, 
we will consider the evolution of the gauge invariant gravitational potential. First consider the 



continuous modes |jT9[. The evolution equation for the gauge invariant gravitational potential $ 
is given by 

„ 6(1-62") , / 2 4(3 + e2'')\ , , 

% - \^%+[p +^- 3_e2. J ^- = 0- (94) 

For small t, $p behaves as T±p{ri)e^^ . For general t, $p behaves as 

Furthermore for small t, this metric perturbation is related to the fiuctuations of the scalar field 
by 

The initial fiuctuations of the scalar field are given in eq.(89). Then, the metric perturbation 
generated during the second infiation is given by 



H 'j:;^ Jo 2v2p V2smh7rp 



+ {e'^^'gl^P. X)T-M + e^P''gl{p, X)Tp{r^))aA +{h.c.)] , (97) 



where 






e^ /i PTl ^2r, 



The variable which has a normalization that relates more directly to the density perturbation 
after reheating is given hy TZ = IQT^GiV'^ /VV)^ . We define the power spectrum of TZ by 

BD{0\n{r,r])n{r',r])\0)BD = J2 dp YpiUr)YpUr')Pn{P,v), (99) 

where a,p\{})BD = 0. Taking the limit rj -^ 0(t -^ oo), we obtain 
Pn{p,X) = Pbd{p)^ 

1 / n — 7 71-1-7 \ " 

,(100) 



12 , I / \M2 - /P-^„ /„ XN„*/„ XA , P + ^ 



\9i{p. W + I^2(P, A)|2 -— ^g,ip, X)g;ip, A) + ^^—glip, X)g2ip, X) 

cosh irp \p + t p — I ^ 

where 

/Si/^y cothvrp 
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The power of the continuous modes in the logarithmic interval p at p ^ 1 is given by 

Here we use the fact the valley bounce gives the initial condition as |0^, — (yj^l ~ M^/ttit, then 
fj,^ = ttitM'^. This quantity should be of the order 10~^° from the observation. This can be 
achieved by taking (M^/M) <C Mp. We show the dependence of A in Bji in Fig. 10. 

The discrete modes can be treated in the same way. $a generated from 6(J)a is given by 

*A = E ^^^^T^{v)YAiUrR, fi), (103) 



where 



e^^ 



Ta = Ta(7i)- 1 + -e'" . (104) 

A ^^'^A + 2 V 3(1 + A) ; ^ ^ 

Taking the limit r^ — > 0, we obtain the power spectrum of TZ 

In some open inflation model, the contribution of the discrete modes gives the strong constraint 
on the model p3[. The supercurvature mode produces very large scale metric perturbations and 



enhances the amplitude of the low multipoles of the CMB anisotropics. No evidence for such 
enhancement in the observed spectrum implies that the contribution of these discrete modes 
must not dominate the contribution of the continuous modes. Furthermore, if the amplitude of 
the supercurvature modes is large, the universe is not open but quasi-open beyond the coherent 
length of the supercurvatuer modes 0|. In our model, however, the last factor in the power 
spectrum P7^(Aj, A) is 0(1) (see Fig. 10), then there is no inconsistency with the observed CMB 
anisotropics and the universe looks like an infinite open universe described by the valley bounce. 

The harmlessness of the supercurvature mode can be deduced from the analysis of the case 
the CD solution describes the tunneling and the thin-wall approximation can be used. In this 
case, the supercritical mode is the wall fluctuation mode given by 50^ = N^{<Pcd)- Here the 
normalization constant is given by A^^ = (/ daa{a)(j)'^j;))^^^^. Within the thin wall approximation, 
this can be evaluated as A^^„ ~ (RSi)^^^^, where Si is the surface tension of the wall and R is the 
radius of the bubble. The surface tension of the wall is estimated by 5*1 ~ ^^(A^)^, where A0 is 
the scale of the variation of the field during tunneling. The curvature perturbation generated from 
this wall fluctuation is given by 7^ ~ {H/(j)cD)S(f)w ~ H/y/RSl. Using R ~ 1/H and rriT ~ H, 
this can be estimated as 7?. ~ H/ /\(j). Thus, the thickness of the barrier the field passes during 
the tunneling determines the amplitude of the curvature perturbation generated from the wall 
fluctuation mode. In the case the valley bounce describes the tunneling, the supercritical mode 
can not be interpreted as the wall fluctuation mode. But the behavior of the supercritical mode 
resembles that of the wall fluctuation mode. Thus we expect this analysis can be applied. Since 
the valley bounce gives A0 ~ M^/rriT, the contribution of the supercritical mode is suppressed. 
Then, the constraint from the discrete mode is not strong in this model. 
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6 Conclusion 



It is difficult to provide the model which solves the horizon problem and at the same time leads 
to the open universe in the context of the usual inflationary scenario. In the one bubble open 
inflationary scenario, the horizon problem is solved by the ffist inflation and the second inflation 
creates the universe with the appropriate Qq. Many works have been done within this framework 
of the scenario and it is recognized this scenario requires additional fine-tuning |0, H . The defect 
comes from the fact that the requirement the curvature around the barrier should be larger than 
the Hubble scale to avoid large fluctuations contradicts to the requirement the curvature of the 
potential should be small to realize the second inflation inside the bubble. Additional constraint 
comes from the fluctuations generated in the decay process, which can be observed and reject 
some models [E3 . 



In this paper we pointed out a possibility of constructing a model without these difficulties by 
reconsidering the tunneling process. If the curvature around the potential is small, the tunneling 
is described by one of a family of the almost saddle-point solutions . This is because the true 
saddle-point solution, that is, the Hawking Moss solution does not satisfy the boundary condition 
for the false vacuum decay. A family of the almost saddle-point solutions generally forms a valley 
line in the functional space. We called the configurations on the valley line valley bounces. To 
identify the valley bounces, we formulated the valley method in the de Sitter spacetime and 
clarified the structure of the valley bounces. In this method the valley bounces can be identified 
using the fact the trajectory they form in the functional space corresponds to the line on which 
the action varies most gently. 

Our assumption is that one of the valley bounces describes the creation of the bubble inside 
of which is our open universe and determines the initial condition in the open universe. Based on 
this assumption, we found that there occurs the second inflation without the large fluctuations 
even if the curvature around the barrier is small compared with the Hubble scale. 

The fluctuations of the tunneling field give rise to the metric perturbation. These can be 
observed in our open universe. The fluctuations around the valley bounce which are orthogonal 
to the gradient of the action are relevant to the observable. We calculated the power spectrum 
of the metric perturbations generated in the second inflation and found these fluctuations can 
be compatible with the observations. In some models of the open inflation, the discrete mode 
of the fluctuations gives strong constraint on the model. We showed this is not the case in 
our model. Hence, using the valley bounce, we can solve the problem which arises in the open 
inflationary scenario besides the usual fine-tuning of the inflationary scenario. Then the one 
bubble open inflation model can be constructed in the simple model with the polynomial form of 
the potential. 

Our conclusion is based on the assumption that an individual valley bounce describes the 
tunneling and gives the initial condition of the tunneling field in the open universe. Although 
this assumption seems to be plausible, we note that further investigations are needed for the 
justification of this assumption. 
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Figure captions 

Fig.l The piece-wise quadratic potential Vt'(0)- Here, we take m^ = 7, rnjp = 2.2 and rj = 0.6M^. 

Fig. 2 The behavior of the CD solution 0(o"), the eigenmode with the negative eigenvalue g-{cr) 
and the scale factor a{a). The potential is taken as in Fig.l. 

Fig. 3 The behavior of the valley bounce. The horizontal coordinate is a. The lower-panel shows 
the behavior of the valley bounce on the trajectory associated with the negative eigenvalue 
and the upper-panel shows the behavior of the valley bounce on the trajectory associated 
with the smallest positive eigenvalue. 



Fig. 4 The action along the valley trajectory. The lower-panel shows the action of the trajectory 
associated with the negative eigenvalue. The horizontal coordinate is A. The upper-panel 
shows the action of the trajectory associated with the smallest positive eigenvalue. The 



horizontal coordinate is the norm of the filed $ = J J daa{a)-^\(f){a) — v^tP- 

Fig. 5 The piece-wise quadratic potential Vt'(0)- We take rri^ = 2, m^ = 0.5 and rj = O.IM^. 

Fig. 6 The action along the valley trajectory associated with smallest positive mode. The poten- 
tial is taken as in Fig. 5. 

Fig. 7 The behavior of the valley bounce in the valley trajectory associated with lowest positive 
mode. The upper-panel shows the behavior of and the lower-panel shows the behavior of 
/. The potential is taken as in Fig.5. 

Fig. 8 Conformal diagram of the de Sitter spacetime. 

Fig. 9 The solution of A. Two solutions are shown. One corresponds to the subcritical mode < 
^sub < 1 and another corresponds to the supercritical mode 1 < A^^p. The corresponding 
valley bounces are shown in fig. 7. 

Fig. 10 The power spectrum of the curvature perturbation around the valley bounce (continuous 
mode). The corresponding valley bounces are shown in fig. 7. 
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